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ABSTRACT 

We defined several functional on the set of all triangulations of the finite 
system of sites in R d achieving global minimum on the Delaunay triangulation 
(DT). We consider a so called "parabolic" functional and prove it attains its 
minimum on DT in all dimensions. As the second example we treat "mean 
radius" functional (mean of circumcircle radii of triangles) for planar triangula- 
tions. As the third example we treat a so called "harmonic" functional. For a 
triangle this functional equals the sum of squares of sides over area. Finally, we 
consider a discrete analog of the Dirichlct functional. DT is optimal for these 
functionals only in dimension two. 

1 Introduction 

Some of the most well-known names in Computational Geometry are those of 
two prominent Russian mathematicians: Georgy F. Voronoi (1868-1908) and 
Boris N. Delaunay (1890-1980). Their considerable contribution to the Number 
Theory and Geometry is well known to experts in these fields. Surprisingly, their 
names (their works remained unread and later re-discovered) became the most 
popular not among "pure" mathematicians, but among the researchers who 
used geometric applications. Such terms as "Voronoi diagram" and "Delaunay 
triangulation" are very important not only for Computational Geometry, but 
also for Geometric Modeling, Image Processing, CAD/CAM, GIS etc. 

The Voronoi diagram is generated by a set of n points S = {p\, in 
R d . The Voronoi diagram is the partition of the R d into n convex cells, the 
Voronoi cells Vt, where each V% contains all points of the R d closer to pi than 
to any other point: 

Vi = {x\Vj ^ i,d(x,pi) < d(x,pj)}, 

where d(x,y) is the Euclidean distance between x and y. 

For generic set of n points S in R d the straight-line dual of the Voronoi 
diagram is triangulation of S 1 , called the Delaunay triangulation and denoted 
by DT(S). The DT(S) is triangulation of the convex hull of S in R d and set of 
vertices of DT{S) is S. 
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Dclaunay triangulation is used in numerous applications. For a plane (d = 2) 
it is usually chosen over other triangulations. Often, it is used in 3D case. A 
logical question may arise: why this triangulation is better than others. Usu- 
ally, the advantages of Delaunay triangulation are rationalized by the max-min 
angle criterion [10, 5]. This criterion requires that the diagonal of every con- 
vex quadrilateral occurring in the triangulation "should be well chosen" [10], in 
the sense that replacement of the chosen diagonal by the alternative one must 
not increase the minimum of the six angles in the two triangles making up the 
quadrilateral. Thus the Delaunay triangulation of a planar point set maximizes 
the minimum angle in any triangle. More specifically, the sequence of triangle 
angles, sorted from sharpest to least sharp, is lexicographically maximized over 
all such sequences constructed from triangulation of S. 

For triangulations of S in 3D and higher dimensions, a sequence of some 
indices corresponding to each simplex can be of different length. Thus, it is 
impossible to compare these sequences. This situation can be circumvented 
if instead of a sequence one considers a functional on a set of triangulations 
S, which provides some real number corresponding to each triangulation. For 
example, instead of minimum angles sequence in max-min criteria, we can con- 
sider a functional which equals to a sum of minimum angles of triangles. This 
functional is minimal on DT(S). 

The purpose of this paper is to present a series of functionals on a set of all 
triangulations of the set 5, which reach their optimum on the Dclaunay trian- 
gulation. These functionals have a clear geometric meaning and demonstrate 
the advantages of the Delaunay triangulation for a plane, and in some cases its 
constraints for higher dimensions. 

2 Optimality of Delaunay triangulations for the 
parabolic functional 

Let t be a triangulation of the set S in R d and let c(Aj) be the center (barycen- 
tcr) of the d-simplex A, i.e. c(Aj) = J^j x ij/(d+ 1), where iy £ R d are vertices 
of simplex A i} j = 0, 1, d. Let 



where |||| is Euclidean norm, and vol(Ai) is the volume (area for d = 2) of the 
simplex A,. 

The functional C2 induces to some order for triangulations of the set S i.e. 
ti > t 2 iff C2(ii) > C2(i 2 )- The value of C2 depends on choice of the origin. If 
we move the origin to x then this order does not change. 

Proposition 1 C2(M) > C2(t 2 ) iffC2(t 1 ,x ) > C2(t 2 ,x Q ). 

The main result for the functional C2 is the following: 



C2(t)=^||c(A i )|| 2 ^(A i ), 
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Theorem 1 The functional C2(t) on triangulations of the set S achieves its 
maximum if and only if t is the Delaunay triangulation. 

Let us consider another functional for triangulations: 



By direct calculation (it is sufficient to verified this formula on a simplex) 
can be proven that 



where CH(S) is convex hull of set S in R d . 

From Theorem 1 and this formula directly follows that: 

Theorem 2 The functional V{t) achieves its minimum if and only if t is the 
Delaunay triangulation. 

A simple proof of these theorems follows from paraboloid construction found by 
Voronoi [11] and rediscovered only in 1979 for a sphere (K.Brown), and later 
also for a paraboloid. 

Let us consider in the space R d+1 a graph of a paraboloid y = \ \x\ | 2 , x 6 R d . 
Let us "lift" the set S from R d into the space R d+1 , x -> (x, \\x\\ 2 ). 

The image of the set 5* creates in R d+1 a set S", whose points lie on the 
paraboloid. The convex hull of the set 5" in R d+1 is divided by vertices which 
were lifted from the boundary of the convex hull of the set S in R d into two 
parts: upper and lower. It turns out that the structure of the polygon of the 
lower part of the convex hull of the set S' when projected on R d gives the 
Delaunay triangulation of the set S. 

Let us consider an arbitrary triangulation t of the set S and "lift" it onto 
the paraboloid in R d+1 i.e. let us build a polyhedral surface in R d+1 connecting 
corresponding vertices on the paraboloid. Note that the functional V equals up 
to a constant to the volume of the solid body below this surface. Thus, the 
minimum V is achieved on the Delaunay triangulation. 

The theorems 1, 2 and the formula relating these functionals show that 
for approximate calculation of the integral: J ||x|| 2 dx with fixed nodes on 



paraboloid, Delaunay triangulation is optimal. This circumstance, of course 
does not prove that the Delaunay triangulation is universally optimal and below 
we will consider other functionals which have a clearer geometrical meaning. 

3 Local circle test 

Let us consider on the plane set of n points S and suppose that any four points 
of that set S do not lie on one circle or line. Let F be functional on the set of 
triangulations of S. 



V(t)=^l + - + ^>ol(A t ). 




CH(S) 



CH(S) 
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For convex quadrilateral ABCD its diagonals give two triangulations: tAc 
and tsD- One of this triangulations is Delaunay triangulation, let it be isn- We 
will say that functional F satisfies Local Circle Test (LCT) if F(tAc) > F(t B D) 
for any convex quadrilateral ABCD. 

Theorem 3 If functional F on the set of triangulations of S satisfies LCT than 
F achieves its minimum on the Delaunay triangulation. 

Actually a proof of this theorem directly follows from [3] , [5] and [2] , when 
using flipping (swapping) algorithm after each flip the functional decreases until 
Delaunay triangulation is reached. 

4 The mean radius of a triangulation 

Let t is triangulation of n points set S on plane. Let us compare each triangle 
Aj of this triangulation with the radius of its circumcircle Ri . Thus every trian- 
gulation t compares with set {i?An --^RAk} of radii of circumcircle of triangles 
Aj E t. The numbers of triangles for any two triangulations of S is equal, so it 
is possible to compare sets of radii for different triangulations. 

In this part we will show that set of radii for Delaunay triangulation is 
"minimal". Here word minimal has concrete definition. In particular, it is 
possible to compare sums of radii: 5Zi?A 4 or power sums: ^R%.,a > 0. It is 
clear that triangulation having minimal set of radiuses is "good", because all of 
its triangles in "medium" the most nearer to the rectilinear triangles. 

Let K is finite set of sequences with k positive numbers: Ri = {Rn, Rik}- 
We call sequence R m as minimal iff f{Rmi) achieves its minimum on K for 
arbitrary increasing function ip. 

Generally speaking, it is possible that for arbitrary set of sequences minimal 
sequence does not exist. One of example, when minimal sequence exists is a 
set of minimum angles of triangles of the triangulation S. In [10] it is practi- 
cally proved that for Delaunay triangulation the sequence of minimum angles 
of triangles is minimal. It turns out that also it is right for sequence of radii. 

Theorem 4 Set of radii of circumcircles of triangles of the Delaunay triangu- 
lation of S is minimal. 

The proof of the theorem follows from Theorem 3 and following lemma: 

Lemma 1 Let ABCD is convex quadrilateral and tsD its Delaunay triangula- 
tion. Then one of the following five relations takes place: 

(1) Rabd < RabcRabd < Radc, Rbcd < Rabc, Rbcd < Radc', 

(2) Rabd < Rabc < Rbcd < Radc', 

(3) Rabd < Radc < Rbcd < Rabc', 

(4) Rbcd < Rabc < Rabd < Radc', 

(5) Rbcd < Radc < Rabd < Rabc, 

where by Rxyz denoted the radius of circumcircle of a triangle XYZ. 
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This lemma shows that for the functional R(t,tp) — X)v 5 (^A i ) satisfies to 
LCT and thus theorem 4 is proved. 

It is an open problem how to generalize this theorem for higher dimensions. 

5 The harmonic index of a triangulation 

The harmonic index of a polygon, a polytope and a triangulation is a pure 
number and prove its name in practice. There are at least two reasons why 
the name "harmonic" takes place here: the first of all this index came from 
so-called "harmonic functions" and the second one is that polygons, polytopes 
and tesselations giving the minimum of harmonic index are more appropriate 
and symmetrical i.e. harmonic. 

For a polygon P its harmonic index is equal the sum of squares of lengths of 
the sides of P divided by area of P i.e. if we denote this index as hrm(P) then 

hrm(P)=J2^/S(P), 

where a\, a n are the lengths of sides of P and S(P) is its area. This index is 
the same for similar polygons. 

We say that an n-gon P as harmonic if and only if hrm(P) achieves its 
minimum for P. It is easy to prove that harmonic triangle is equiangular. The 
same result hold for arbitrary n. 

Proposition 2 A polygon P is harmonic if and only if P is a regular polygon. 

The proof of this proposition very close to proof of minimal circle properties 
(more specifically, isoperimetric inequality for polygons) in book [1]. 

It is easy to calculate hrm for a harmonic n-gon H, hrm(H) = 4tan(7r/n), 
i.e. for arbitrary n— gon P, hrm(P) > 4tan(-7r/n). 

For a planar triangulation t of sites S let denote by hrm(t) (the harmonic in- 
dex of the triangulation t) the sum of hrm of its triangles: hrm(t) — hrm(Ai) 

Our main result for the harmonic index of planar triangulations is the following: 

Theorem 5 The harmonic index hrm(t) of a triangulation t of S achieves its 
minimum if and only if t is the Delaunay triangulation of S. 

In other words, the set of harmonic indices of DT triangles is minimal. This 
theorem follows from theorem 3 and fact that hrm(DT) satisfies to LCT. 

The harmonic functional for triangle gives its minimum if this triangle is 
equiangular. Usually, a triangulation is regarded as "good" for different pur- 
poses if its triangles are nearly equiangular. The harmonic index of triangulation 
t achieves its absolute minimum if t is part of a regular triangular lattice. In 
some sense, theorem 5 shows that the Delaunay triangulation among all trian- 
gulation of S as so close to equiangular triangulation as possible. 

We have proposed the generalization of the harmonic index for polygons. 
Let for polygon P: hrm(P,k) — ^af fc /S'(P) fc , where k is some real number. 
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If k > 1/2 then hrm{P, k) achieves its minimum for rectilinear polygon. In the 
case k = 1/2 this result gives classic isoperimetric inequality for polygons [1]. 

It is easy to give an example of convex quadrilateral P when its triangulation 
t giving minimum of functional hrm(t, k),k ^ 1 is not Dclaunay triangulation. 
In this connection open problem appears: "For fixed k > 1/2 find efficient 
algorithm for construction of such triangulation t of S, that functional hrm(t, k) 
achieves its minimum." We know the answer only for the case k = 1. It is any 
of algorithms of construction of Delaunay triangulation. 

Let us consider the harmonic index in d dimensions. It is clear that right 
extension of hrm for polytope P is the following: 

hrm(P)=J2F?/Vol d - 1 , 

where Fi are the volumes of faces of P and Vol is the volume of P. 

As in 2D case a harmonic tetrahedron (simplex) is regular. It seems to us 
that all Platonic solids are harmonic. 

It is easy to construct d+2, d > 2 vertices polytope that its Delaunay trian- 
gulation does not gives minimum of hrm. Open problem is to study "harmonic" 
triangulations and to find an efficient algorithm for its construction if d > 2. 

6 The Dirichlet functional on triangulations 

Let S — (xi, ...,x n ) be a set of points in R d , each associated with a real number 
Hi. Denote by Y the set of these numbers, i.e. Y = (y\, ...,y n ). There are a 
lot of different problems in GIS, Geology Topography CAD/CAM etc., where 
need construct a surface in R d+1 corresponding to this set of data. The main 
problem is following: to find a function y = f(x), such that f{xi) — yi. One 
of the oldest and most famous methods is modeling by triangulation. If we 
have some triangulation of S then for a set of data Y there is only one method 
to construct a piecewise linear function (polyhedra) on this triangulation. The 
question is what triangulation is "good" for modeling. Usually, for this purposes 
to use the Delaunay triangulation. 

One of the criterions that triangulation is suitable for the modeling is that 
fact that area (volume if d > 2) of the polyhedral surface y = f{x)) in R d+1 is 
minimum. Let us introduce the following functional: 

SV(t,Y) = Y / vol(A t (Y)) 7 

i 

where (Ai(Y) denotes simplex on the surface y = f(x) in R d+1 over the (A*. 

It turns out that on the plane for " small" yi triangulation, giving minimum 
for this functional does not depend from Y and it is DT. 

Theorem 6 For every set of points S in the plane there is real number e > 
such that if for all i : yi < e, then functional SV(t, Y) achieves its minimum iff 
t is Delaunay triangulation. 
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We can give an example of 6 points on the plane and some set Y when the 
minimum of the functional SV not reaches on the Delaunay triangulation. 

In this connection for the arbitrary set Y open problem appears: to find an 
algorithm for constructing triangulation that gives minimum for functional SV 

The other minimum criteria is a discrete analog of the Dirichlet functional: 
/ grad 2 f(x)dx. For interval (d = 1), a spline deg — 2k — 1 is function y = f(x) 
such that f(xi) = yi and 

b 

J [f {k \x)] 2 dx = min. 

a 

It is clear that discrete analog for k = 1, d > 1 and piecewise linear function 
/ is: 

DF(t,Y)= J \\gradf(x)\\ 2 dx = ]T ^ - vol(CH(S)). 

CH{S) 1 

Triangulation where this functional reaches minimum is possible to call dis- 
crete spline triangulation (DST). For the plane DST does not depend from Y 
and it is DT. 

Theorem 7 (S.Rippa [9]) DF(t,Y) achieves its minimum ifft is the Delaunay 
triangulation. 

The proof of the theorems 6 and 7 follows from theorem 3 and fact that this 
functionals satisfies to LCT. It is possible to prove it directly. It also follows 
from some general result that is given bellow. 

Let S be a set of d + 2 points x\, Xd+2 in R d . C.Lawson [6] proved that 
there are at most two distinct triangulations of convex hull of S. Suppose S 
admits two triangulations t\ and t 2l and Y = (yi, yd+2) is a set numbers 
corresponding to x\, Xd+2 as above. Let B(Y,S) = DS(t\,Y) — DS(t2,Y). 
Note B(Y, S) is a quadratic form depending on Y. 

Theorem 8 The optimal (DST) triangulation for n = d + 2 does not depend 
on Y 

This theorem follows from the fact that B{Y) = const(S)L 2 (Y), where L(Y) 
is some linear form on Y. 

For d > 2 is easy to give example when triangulation is DST, but not DT. 
Open problem is: To study discrete spline triangulation for d > 2. Can it 
depends on Y if n > d + 2 ? 

7 CONCLUSION 

Functionals given in this paper could be useful in analysis of algorithms and con- 
struction of optimal triangulations. For example, the proof that after sequences 
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of " local transformation" in [4] we do not get initial triangulation easily follows 
from consideration of the functional V(t). The point is that this functional 
decreases after each local transformation, used in the algorithm. In [4] quite 
complicated proof is given for it. 

Analysis of functional " the mean radius" , hrm and DF demonstrates that 
for d > 2 the Delaunay triangulation is not optimal. It makes doubts that 
for applications in higher dimensions the Delaunay triangulation can be used 
everywhere. The problem of finding "good" triangulations for these functionals 
in higher dimensions is open and so more detailed consideration is necessary. 
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